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Abstract 

This work adresses the question of the scaling behaviour of poly- 
electrolytes in solution for a realistic prototype: We show results of 
a combined experimental (light scattering) and theoretical (computer 
simulations) investigation of structural properties of poly (acrylic acid) 
(PAA). Experimentally, we determined the molecular weight {Mw) 
and the hydrodynamic radius (Rh) by static light scattering for six 
different PAA samples in aqueous NaCl-containing solution (0.1 — 1 
mol/L) of polydispersity Dp between 1.5 and 1.8. On the computa- 
tional side, three different variants of a newly developed mesoscopic 
force field for PAA were employed to determine Rh for monodisperse 
systems of the same Mw as in the experiments. The force field effec- 
tively incorporates atomistic information and one coarse-grained bead 
corresponds to one PAA monomer. We find that Rh matches with the 
experimental data for all investigated samples. The effective scaling 
exponent for Rh is found to be around 0.55, which is well below its 
asymptotic value for good solvents. Additionally, data for the radius 
of gyration (Rg) are presented. 
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1 Introduction 



Poly (acrylic acid) (PAA) is a water-soluble polyelectrolyte (PE) . It is impor- 
tant not only in industrial applications, e.g. flocculants or superabsorbers.0 
Because of its relatively simple chemical repeat unit, it is also a prototype 
PE model for scientific investigations. Surprisingly, there seems to be no (ac- 
curate) data on some solution properties of PAA chains. Quantities which 
characterize the solution behaviour of isolated polymer chains are the radius 
of gyration Rq and the hydrodynamic radius Rh- It is the main purpose of 
this letter to provide reliable values for PAA of different chainlengths. 

Two approaches are being used. Experimentally, we determine the size of 
a polymer coil in solution by dynamic {Rh) or static {Rg) light scattering. 
The measurements were performed using narrow fractions of radically poly- 
merized PAA in order to, on the one hand, minimize effects of polydispersity 
and to, on the other hand, avoid possible problems of PAA aggregation due 
to hydrophobic initiators used in anionic polymerization.! These data are 
augmented by results from computer simulations. Simulations have the prin- 
cipal advantage that macroscopic observations can be understood in terms of 
a microscopic model. We will use them here to investigate the scaling of the 
polymer also at length scales shorter than the overall coil size {Rq or Rh) 
which is accessible by the light-scattering experiment. This is the second aim 
of this work. 

The experimentally relevant molecular weights are far beyond what can be 
simulated with an atomistic model. They are, however, accessible with suit- 
ably simplified or "coarse-grained" (CG) models. Care has to be taken that 
the CG model retains sufficient information about the chemical nature of the 
polymer. A generic bead-spring model might be good enough for theoreti- 
cal scaling relationsjl^ but will not give realistic absolute values for Rh of 
specific polymers like PAA. A number of systematic coarse-graining proce- 
dures have been described,IHl§ but none has been tried for specific PEs in 
solution. We have recently developed an automatic coarse-graining method 
which works by a two-step process: Firstly, an atomistic simulation of a PAA 
oligomer in water is performed. Secondly, the CG model is parametrized to 
reproduce the PAA structure.0 To validate this method and to try out 
technical variants on it is the third objective of the letter. 
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2 Methods 



2.1 Experimental 

2.1.1 Static and Dynamic Light Scattering 

For simultaneous static and dynamic light scattering measurements a com- 
mercial instrument (ALV-5000) with a krypton ion laser operating at a wave- 
length of 647 nm and an Avalanche Diode (EG&G) as detector was used.lil 
Static measurements were performed at scattering angles of 30° — 150° in 10° 
steps in the concentration range of 0.5—5 g/L. For the high molecular weights 
the highest concentration was only 3.5 g/L in order to remain in the dilute 
regime. Weight average molar masses and radii of gyration Rq were 
obtained by Zimm extrapolation using the Rayleigh ratio RR = 1.270 cm~^. 
The refractive index increments dn/dc are listed in Table ^, as measured at 
A = 663 nm using a scanning Michelson interferometer .Ej 

Dynamic measurements were performed over an angle range from 60° — 150° 
in 30° steps. The intensity autocorrelation function of the scattered intensity 
was converted to the scattered electric field and analyzed using the program 
CONTIN by S.Provencher.0 Hydrodynamic radii Rh were obtained via 
the Stokes-Einstein relation, where the apparent diffusion coefficient D = 
T /(f was calculated from the inverse relaxation time F and the absolute 
value q of the scattering vector .El Viscosity corrections due to the NaCl were 
considered to be small and not taken into account. The hydrodynamic radius 
Rh was measured repeatedly. In Figure |^ the average values and error bars 
representing one standard deviation are plotted. 

2.1.2 Sample Preparation 

The PAA samples (Polymer Standard Service) with molar masses of My/ = 
18000 — 296600 g/mol were synthesized by radical polymerization and frac- 
tionated with polydispersities Dp = Mw/Mn between 1.5 and 1.8 (with Mn 
being the number average molecular weight). All samples were dissolved 
in deionized water (Millipore) with IM NaCl. The pH was determined for 
two samples {M\y = 81800 and 296600) at two different concentrations (1 
g/L and 5 g/L) and found to lie between 6 — 8. After further purification, 
we performed additional measurements by dialysis using a regenerated cel- 
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lulose membran (MWCO2000). The dialysis process lasted 40 hours and the 
deionized water (Millipore) was replaced four times. Afterwards, the sam- 
ples were freeze-dried for 5 hours and finally dried under vaccuum for 60 
hours at T = 303 K, in order to avoid fluctuations of the PAA concentra- 
tion by adsorbed water. The refractive index increments for the dialyzed and 
dried samples remained unchanged within the statistical errors. All solutions 
were filtered through a 0.22 fim Millex-GS filter (Millipore) to remove dust 
particles. 



2.2 Computational 

2.2.1 The Coarse-grained Force Field for PAA 

Coarse-grained potential energy functions for polymers have to incorporate 
not only energetic, i.e. local aspects of the underlying microscopic model. 
They also have to account for entropic contributions from the neglected con- 
formational degrees of freedom of the chain. Therefore, we utilized structural 
information obtained from fully atomistic simulations as target functions to 
construct our CG force field. In particular, we used the intra- and inter- 
chain radial distribution functions which, being distributions derived from 
an ensemble, contain the desired entropic information. This is a so-called 
inverse problem: find an interparticle potential which reproduces a given 
radial distribution function (RDF) or set of RDFs. For the fitting proce- 
dure, we applied an automatic optimization algorithm (simplex) which was 
originally implemented for the development of atomistic force fieldii and 
tested for liquids.ii Our force field is based on atomistic simulation data 
by Biermann et al@ They studied one fully deprotonated, atactic oligomer 
of 23 monomers with 23 Na"*" counterions in 3684 water molecules (simple 
point charge (SPC) modelil) at ambient conditions. The system represents a 
diluted (~ 2.3wt%) PAA solution. From this simulation, structural informa- 
tion like the distributions of bond angles or RDFs between monomers were 
extracted. We mapped this system to the mesoscale by replacing each repeat 
unit (i.e. each monomer) by one bead, either at the monomers center of mass 
or at the backbone carbon bearing the carboxyl group. The PAA mesoscale 
force field includes bonded and non-bonded interactions. To both contribute 
several terms.i^ They were parametrised by systematically varying the in- 
teractions until the structure of the atomistic model was reproduced. This 
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also allowed us to omit all explicit water molecules and counter ions. Their 
effect on the PAA chain conformation is, however, implicitly present in the 
model. In effect, a system of roughly 10000 Atoms could thus be reduced to 
a system which consists of only 23 "superatoms". 

For the sake of gaining experience with CG force fields, we have coarse- 
grained PAA in three different ways. Potential / is pieced together from 
various analytic functions. It contains attractive as well as repulsive inter- 
molecular interactions. Coulomb forces were only implicitly taken into ac- 
count (they were present in the parent atomistic simulation). That is because 
the Debye length (characterizing the relative importance of the electrostatic 
energy compared to thermal energy) of the system is close to that of water, 
so that electrostatic interactions can be treated as being short-ranged. For 
details, see ref 17. Potential // is similar in construction, but uses another lo- 
cation of the CG bead: While Potential I placed it at the center of mass of a 
monomer. Potential // employed the backbone carbon parent to the C00~- 
group.0 As example for the achieved convergence, the RDF (first neighbours 
excluded) is shown in Figure |T]a. Potential III uses again the center-of-mass 
of the PAA monomer, but was optimized with a fully numerical instead of 
a piecewise analytical potential energy function for the non-bonded interac- 
tions.il It allows us to apply a new self-consistent optimization algorithm 
which was recently adapted by usJ^l The optimization then yields an almost 
perfect match of the RDF. Target RDFs (exluding first and second neigh- 
bours, as defined in Ref.0) and the results of the optimization are shown in 
Figure ^d. The parameters of force fields / and // are presented in Table |1|. 
The non-bonded part of Potential III exists in numerical form onlyS (not 
shown here). 

2.2.2 Technical Simulation Details 

Both Brownian Dynamics (BD) and Monte Carlo (MC) programs were used 
to carry out the PAA simulations. We simulated at T = 333.15 K, which 
corresponds to the temperature of the atomistic simulations. The Langevin 
equations of motion were integrated by the velocity Verlet algorithm with 
a time step At = O.Olr,!! and a friction constant 7 = 0.5r"^xl Pivot MC 
calculations were necessary to simulate systems of more than 1000 monomers 
as the equilibration is much faster. We carried out 10^ accepted equilibration 
moves before a production run of 10^ accepted moves started. 
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3 Results and Discussion 



The hydrodynamic radii of six different PAA samples witfi molecular weights 
in the range from 18100 to 296600 g/mol were measured (see Table |]). For 
four samples, the molar masses Mw and the radii of gyration Rq were mea- 
sured in IM NaCl solution. (For the samples of lower molecular weight an 
accurate determination of Rq was not possible, so they are reported as speci- 
fied by the supplier.) Additionally, the ratio Rq/Rh and the refractive index 
increment dn/dc are shown in Table ^ In all PAA solutions, the salt con- 
tent was so high that the electrostatic interactions were effectively screened. 
Therefore, Rq could be obtained from the static scattering experiment by 
Zimm extrapolation to zero angle and zero polymer concentration Cp like 
for uncharged polymers. The dynamic measurements in IM NaCl solution 
yielded single-exponential autocorrelation functions, since we performed our 
measurements in the "ordinary regime" where c™/c™ < 1, with c™ and c™ 
being the molar concentration of monomer units and of salt, respectively.^^ 
No polyelectrolyte slow mode due to a cooperative diffusion behavior was 
observed. Furthermore, the diffusion coefficients do not show any significant 
angular dependence. 

For two samples with molecular weigths of = 119800 g/mol and M\y = 
204000 g/mol, we varied the NaCl concentration between 0.1 — 1 mol/L at a 
fixed polymer concentration of 1 g/L, as shown in Figure |^a. No dependence 
of the diffusion coefficient or Rh is found within this range. Figure ^ shows 
a weak decrease of Rh with polymer concentration in the range of — 5 g/L 
for two molecular weights at fixed NaCl concentration of 1 M. For fur- 
ther comparisons we use the values extrapolated to zero concentration. Thus 
the static as well as dynamic light scattering experiments represent mea- 
surements of single polylelectrolyte molecules. The hydrodynamic radii are 
plotted in Figure ^ as a function of molecular mass. We find a scaling ex- 
ponent of z/ = 0.56 ± 0.02. This scaling exponent compares well with CG 
simulation results (below). 

In Table ^, the corresponding theoretical data is presented for 13 different 
chain lengths (being single chains, the simulated polymers are monodisperse). 
They include all the experimental ones when converting the mean molar 
weights Mw into the degree of polymerization. We plotted them in Figure |^a, 
too. Over the whole range of measured samples, the coincidence is excellent. 
This is especially remarkable if one considers the fact that the simulation 
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model was developed with a PAA chain consisting of only 23 monomers. 

In order to check the transferability of parameters between oligomer chains 
of different lengths, we executed new atomistic simulations of shorter chains 
with the same force field as used by Biermann et al. in their work.@ The 
PAA was dissolved in 4000 water molecules (concentration Cp ~ 1 g/L) and 
the total simulation time was 6 ns. The concentration of sodium counterions 
was determined to be around 0.44 mol/L for all atomistic simulations, which 
lies in the middle of the experimentally tested ones. A "downward" trans- 
ferability could be obtained for both samples. The CG model developed for 
a 23-monomer chain, faithfully reproduces also the Rh and Rq of the atom- 
istic 8- and 12-mers, cf. Table |^. It is also interesting to investigate the 
stability of Rh for the various simulation models /, // and ///. As shown in 
Table ^d, the coincidence is good for all chain lengths. This means that our 
CG mapping is not unique. Structural properties can be determined from 
any of the above model variants. That might be of use, if one is interested in 
properties pertaining to the center-of-mass or of the CH-backbone carbon. 

In the following, we focus our attention on the original force field I. The sim- 
ulated Rh for chains longer than 100 monomers was fitted by a power law 
Rh oc M^. The resulting scaling exponent v was 0.55 ± 0.02. If we restrict 
the fit of the experimental data to the simulated chain lengths, we obtain 
the same exponent. Experiments and simulation yield, hence, very similar 
results. It can be concluded that the salt concentration (used in experiments) 
or sodium counterion concentration (in the parent atomistic simulations) ef- 
fectively screen the charges of the PAA chain. Thus, the Coulomb potential 
of PAA is really reduced to an effective short-ranged interaction. For small 
distances, though, the charges stiffen the chain. This behaviour can be un- 
derstood in the framework of scaling theory. For a Gaussian chain (model of 
a linear fuUy-fiexible polymer) in a good solvent, self-avoiding walk (SAW) 
statistics applies. It vields a scaling exponent of = 0.588 in the limit of 
infinitely long chains.a This scaling factor is universal, i.e. it is valid for ev- 
ery size-characterizing function, e.g. Rq or Rh- The structur factor S{q), 
however, is an even better function to investigate. That is, because it scales 
as S{q) cx q'^^" 1 which we can try to fit for any length scale up to the total 
size of the polymer. So, local deviations from the global behaviour, as mea- 
sured by the experiment, can be observed. Figure ^ shows the structure 
factor of the coarse-grained BD simulation of a PAA 460-mer. Apart from 
noise in the region q > 10, three different regimes can be distinguished: On 
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very short distances (nearest neighbours), the slope equals —1.65 = — 
i.e. corresponding to a "local" scaling exponent of u = 0.61. This is remark- 
ably close to, but larger than the value for a SAW (0.588). Compared to 
a Gaussian chain, this is due to the stiffening bond angle potential of the 
chain (in the CG picture). Next comes the region in which torsions and non- 
bonded interactions become important. They introduce even more stiffness, 
resulting in a large exponent of u = 0.82. Finally, the long-range behaviour 
sets in (z/ = 0.56), almost yielding the obtained total scaling exponent for 
Rh- So we can conclude that PA A globally behaves like a fully- flexible chain 
without charges. These results are very similar for any other large-enough 
PAA sample. We note also that the original target RDF of the 23-mer is 
reproduced by any 23- monomer sub-strand of the 460-mer (data not shown). 
For Rg, scaling similar to that of Rh is, as expected, obtained in the limit 
of sufficiently long chains (> 100 monomers). In Figure this can be seen 
by examining the ratio Rg/Rh'- a plateau region is reached for long chains, 
with a value of around 1.5. This is the value predicted for Gaussian chains. 
A closer look, however, reveals a monotonically increasing ratio (cf. Table 
exceeding the value of 1.5. The reason for this behaviour are the corrections 
to scaling due to the good solvent statistics of the chain, which have been 
known for a long time.o For Rh, this motivated the empirical indroduction 
of an effective exponent u^ff ~ 0.55,0 which exactly equals our findings. A 
detailed analysis of these corrections, whose asymptotic value is found to be 
Rg/Rh ~ 1.61, will be presented elsewhere.0 
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Table 1: Parameters for our various coarse-grained force fields. For details, 
see refs 17 and 27. 



(a) bond angle parameters 

Peak # position standard deviation peak height 

[°] [°] [height of peak 2] 



center of mass 
force fields /,/// 

1 88.0 6.8 1.72 

2 116.7 7.4 1.00 



backbone carbon 
force field // 


1 
2 


127.4 
156.4 




8.5 
8.9 




1.15 
1.00 




(b) non-bonded parameters 


















0-1 


£i 0-2 






£3 




64 


original force field 
(center of mass) / 


0.496 


11.30 0.559 


0.35 


0.775 


0.49 


1.3 


0.05 


backbone carbon 
force field II 


0.496 


11.30 0.565 


1.20 


0.802 


0.70 


1.3 


0.05 



Table 2: Experimental hydro dynamic radii Rh and radii of gyration Rq 
obtained by dynamic and static light scattering experiments for poly (acrylic 
acid) of different molecular weights Mw- 



Mw 


g/ mol] 


Rc [nm 


Rh [nm 


Rg/Rh 


dn/dc 


18100^ (193-mer) 
36900^ (393-mer) 
81800^ (870-mer) 
119800^ (1287-mer) 
204000^ (2170-mer) 
296600=^ (3155-mer) 


12.9 
14.9 
22.7 
23.8 


3.7 

5.3 

8.9 

10.4 

15.0 

16.6 


1.45 
1.44 
1.51 
1.46 


0.1360 
0.1452 
0.1472 
0.1310 



^: Mw as specified by the supplier 

2; Mw determined by own static light scattering 
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Table 3: Radii of gyration Rq and hydrodynamic radii Rh for different- 
size poly (acrylic acid) as obtained by Brownian Dynamics (marked '1') and 
Monte-Carlo (marked '2') simulations. In the coarse-grained picture, poten- 
tials / and /// are different versions of center-of-mass mapped force field, 
whereas potential // is centered at the backbone carbon of the chain. 



Chainlength 


Rg [nm] 


Rh [nm] 


Rg/ Rh 


(a) Atomistic Model 


8-mer 


0.55 


0.87 


0.63 


12-mer 


0.77 


0.91 


0.85 


23-mer 


1.28 


1.19 


1.08 


(b) Coarse-Grained Model 


Potential / 








8-mer^'^ 


U.O ( 


u.oy 


n P.A 

U.D4 


12-mer^ 


0.78 


0.94 


0.83 


23-meri'2 


1.25 


1.18 


1.06 


46-mer^ 


2.0 


1.6 


1.25 


92-mer^ 


3.1 


2.3 


1.35 


193-mer^ 


4.8 


3.4 


1.41 


393-mer^ 


6.8 


4.8 


1.42 


460-mer^ 


7.7 


5.3 


1.45 


460-mer^ 


7.9 


5.5 


1.44 


852-mer^ 


10.8 


7.4 


1.46 


855-mer^ 


11.3 


7.7 


1.47 


1287-meri 


13.5 


9.1 


1.48 


2067-mer2 


18.8 


12.4 


1.52 


3155-mer^ 


24.1 


15.7 


1.54 


Potential // 








23-mer^ 


1.22 


1.03 


1.18 


393-mer^ 


6.9 


4.8 


1.43 


852-mer^ 


9.4 


6.6 


1.42 


1287-mer^ 


14.3 


9.0 


1.59 


Potential III 








23-mer^ 


1.26 


1.18 


1.07 


460-mer^ 


8.6 


5.9 


1.46 


852-mer^ 


12.2 


8.3 


1.47 


1287-mer^ 


14 15.5 


10.3 


1.51 



60 




Radius r [nm] 
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